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Abstract 

In this study we have analytically obtained relaxation function in terms of rotational correlation 
functions based on Brownian motion for complex disordered systems in a stochastic framework. We 
found out that rotational relaxation function has a fractional form for complex disordered systems, 
which indicates relaxation has non-exponential character obeys to Kohlrausch- William- Watts law, 
following the Mittag-Lemer decay. 



PACS numbers: 05.40.-a, 02.50.-r, 02.30.-f, 76.20.+q 



Relaxation for ordered systems is given by Maxwell-Debye law as 

= $ exp(-t/r), *>0. (1) 



Whereas, re 



axation in many complex disordered systems such as metallic glasses, spin 



glass alloys ^, ^ , 0, |(| , ferroelectric crystals [?J , dielectrics deviates from the classical 
exponential Maxwell-Debye pattern and is often described in terms of Kohlrausch- William- 
Watts (KWW) (i.e. stretched exponential) law 0,0], 

$(t)=$ exp(-i/rr (2) 

for < a < 1, or by an asymptotic power-law 

$(t) = $ (l + t/rp (3) 

with n > 0. The relaxation functions in Eqs.(l) and Eq.(2) are commonly written in terms of 
the correlation functions which correspond to decay of the fluctuation of a physical quantity 
such as magnetization in magnetic materials or polarization in the dielectric materials. 

Relaxation function has been derived using rotational relaxation method for some sys- 
tem. This method has been firstly used by Debye in the context of dielectric relaxation of 
polar molecules Q . Debye theory is based on the Smoluchowski equation for the non-inertial 
rotational diffusion of the molecules. In his work on dielectric relaxation of an assembly of 
non interacting dipolar molecules Debye considered two models of the process, namely (a) 
an assembly of fixed axis rotators each having a permanent dipole moment /x and subjected 
to Brownian motion torques having their origin in the background or heat bath and (b) 



the same assembly however the restriction to fixed axis rotation is removed. T 
in both instances are equivalently the same, if inertial effects are disregarded 



re results 
The 



3. 



same picture has been applied to the rotational motion of the magnetization vector of a 



super-paramagnetic particle 



11 L the polarization vector of a polar molecule in a dielectric 



j!0|. and heavy molecules [12. 13] in the liquid and gases. However. Debye theory cannot 
explain the experimental data on dielectric relaxation of complex systems, since the inter- 
actions between dipoles are ignored. Indeed, the relaxation process in disordered systems is 
characterized by the temporal nonlocal behavior arising from the energetic disorder which 
produces obstacles or traps which delay the motion of the particles and introduce memory 
effects into the motion. Therefore, an important task in relaxation of complex systems as 



well as dielectric relaxation is to extend the Debye theory of relaxation to fractional dy- 



namics, so that empirica 
y, Williams and Watts 

n r 

was given in Refs. \14L |15 



decay functions, e.g., the stretched exponential of Kohlrausch 
£ I, may be justified. Such a generalization of the Debye theory 
, 1^ . We must remark that other generalizations of the Debye 



model in the context of the fractional dynamics have been discussed in recent striking works 



17 



18 



y 



20 



2l| as well. 

Our main aim, in this study, is obtain relaxation function Eq.(2) for some complex dis- 
ordered systems in terms of rotational correlation functions based on rotational Brownian 
motion in a stochastic framework. 

The most common calculation in which the picture of rotational Brownian motion finds 
relevance is that of the rotational correlation function which can be measured by 

Infrared (IR) and Raman (R) spectroscopies as well as the neutron and ultrasonic scattering 
techniques |22j. Correlation function <3>f (t) measures the correlation in time between the 
direction of the unit vectors (defining the molecular symmetry axis) u(0) and u(t). There- 
fore, $f (t) as can be obtained from the IR and Raman spectroscopies may be expressed in 
the compact form, 

d>f (t) = (fl(u(0).u(t))> (4) 

where u is the unit vector along the symmetry axis of the molecule, and Pi is the Legendre 
polynomial of order I. Argument of Pi is expressed as u (0) .u (t) = cos 7 (t), where 7 is angle 
between two different point at the spherical coordinate system. These points are given by 
{9o, 4>o) an d (9, 0) which denote two different directions by separated by an angle 7. There 
angles satisfy the trigonometric identity, 

cos 7 = cos 6*o cos 9 + sin 9q sin 9 cos (0o — <P) (5) 
The addition theorem asserts that 

4-7T ^ 

Pi (cos 7) = ^ Yl Y i m (*o, <M YT (0, 4>) ■ ( 6 ) 

m=—l 

If we inserted Eq.(6) into Eq.(4) rotational correlation function is expressed in general form 
using spherical harmonics addition theorem as 

4-7T 

$1 (t) = ^7 £ (Yr (0o, 0o) Yr (9, <P)) . (7) 

m=— I 
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Such a theoretical approach allow us to calculate for arbitrary I. Average in Eq.(7) is 
calculated using the probability density functions of Brownian particle which refers to end 
of point polarization vector or a real particle. The Brownian motion in disordered space in 



the presence of an externa 



(FFPE) 
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24, 



25 



26 



2J: 



field F(x) = —V'(x) leads to Fractional Fokker-Planck Equation 



^-W (x , | x,t) = D]- a L FP W (x , | x,t) 



(8) 



jfp 



d V (x) T ^ d 2 1 
+ K a — 



(9) 



dx mrj a dx 2 

This equation then characterize sub-diffusion process. The FFPE are closely related general- 



ized Le vy-t ype statistics 



models 



29, 3( 



31 



28 [and can be derived from continuous time random walk (CTRW) 



34|, or from a Langevin equation [35(. In Eq.(8), W (x ,0 \ x,t) 



imply conditional probability for Brownian motion, m denotes mass of the particle, K a the 
diffusion constants associated with the transport process, and the friction coefficient r] a is a 
measure for interaction of the particle with its environment. K a is a generalization of the 



Einstein-Stokes-Smoluchowski relation 



24, 



25 



26 



271 ] which holds for the generalized 



coefficient r] a , which is defined as K a = kBT/mr) a where kp is the Boltzman constant, and 
T is the temperature. In theJEq.(8) operator ofj"" is known fractional Riemann-Liouville 
integro-differential operator 



The fractional Riemann-Liouville operator Q D\ a 

i 

1 d 



Q D t a featuring is defined through 



r)l-a 



W (flo,0o,O I 9,<f>,t) 



T (a) dt 



dt 



,W(0 o ,</>o,O|0,M 



(t - 1') 



l-a 



(10) 



The fractional integro-differentiation operator qD\ a contains a convolution integral with a 



slowly decaying power-law Kernel M (t) 



r(a) 



ensures the non-Markovian nature of the 



sub-diffusion process defined by the fractional diffusion process. Its fundamental property 
is the fractional integro-differentiation of a power, 



p,i-arp _ r (1 + p) +p+a _i 
T (p + a) 

In fact, it can be shown that more general relation 

r(i + <z) 



(ii) 



D p A q 



T(l + q-p) 



t g- P 



(12) 



for any real p, q. Thus, the fractional derivative of a constant, 

1 



-r 



q > 



(13) 



r(l-g) 

reproduces an inverse power-law. The special cases of integer order integro-differentiation of 
a constant, = 0, are included through the poles of the Gamma function for q = 1, 2, 3, .... 

In the case of V (x) = which means that there is no an external field, the one- 
dimensional FFPE can be reduced to diffusive type equation; 



-W (z , | x,t) = D]- a K a —W (x , | x,t) . 



37 



38 



(14) 

which is a particular form of 



This equation is called fractional diffusion equation 
the FFPE, which can be represent on the spherical coordinates as a function of the angles 
9 and 6 as 



W (9 , 0o, | 0,0, t) = D]~ a d a W 2 W (0 O , O , | 9, <j>, t) 



(15) 



where d a is referred to as the rotational diffusion constant that is related to translational 
diffusion constant K a and the radius a by d a = K a /a 2 (a = unit) 



sin 2 9 



Smd d9[ Smd d9 



+ 



d 



2 1 



(16) 



The standard method of solution of Eq.(15) is the separation variables j3, 2d, 13 Q- 
271 ] . If we consider separation ansatz, as W = T.Q, where T and Q temporal and spatial 
components of the conditional probability function respectively, we will obtain two eigen- 
equations as 

. rr i f\ I i \ 

\ 2 d aQ D l t - a T{Q\t) (17) 



dT (0 | t) 



dt 



V 2 Q (0 , 0o \9,<p) = -\ 2 Q (9 , O | 9, 0) (18) 

The temporal eigen-equation Eq.(17) is but the fractional relaxation equation, the solution 
of which is given in terms of the Mittag-Lefner function 

-dj {I + 1) t a ] j 



T(Q\t) = E a [-d a l (l + l)t a ] = J2 



3=0 



r (1 + aj) 



(19) 



As can be seen from the series expansion, the exponential form can be recovered in the 
Brownian limit a = 1 



E x [dxl ( I + 1) t 1 ] = exp [dxl ( I + 1) t] . 
5 



(20) 



This result indicates that for a = 1 Mittag-Leffler relaxation modes leads to Maxwell-Debye 
pattern, on the other hand, for interval < a < 1 Mittag-Leffler function indicates non- 
Markovian Brownian motion. 

On the other hand, a convenient way of expressing solution of the spatial eigen-equation 
Eq.(18) is by means of spherical harmonics 

oo I . 

Q (0 , 0o I M) = EE 27TT^ m ^ ^ Yr ^ ■ (21) 

1=0 m=-l 

As a result, if solutions Eqs.(19) and (21) are combined, conditional probability function 
of Brownian particle is obtained in terms of Mittag-Leffler function as 

oo I . 

w (e , 0o, o | M,t) = EE 2zTI y,m ( ^ ' 0o) Yr {6, 0) Ea [ ~ dJ {l + 1} r] (22) 

1=0 m=-l 

Eq.(22) can be written as a compact in terms of the Euler angles n o (#o,0o) and f2(6>,0), 
hence 

oo / . 

w (no, o | n,t) = J2 zrri Y ™ W Y ™ ^ E « ^ dJ + ^ ( 23 ) 

i=0 m=-l 

Also, Eqs.(22) and (23) should be satisfied 

w(n ,o | n,o) = 5(n -n) (24) 

for i = 0. 

It is easy give the correlation function Eq.(7) what is by now a familiar interpretation as 

A ^ P 

m=— I 

In above equation w (Qq) is a priori probability that the initial orientation is given by Qq, 
while (Qo, | f2, t) is the conditional probability that the final orientation is determined Q. 
Assuming that the reorientations of the spin (or molecular) symmetry axis may be modelled 
as an isotropic rotational Brownian motion, we may write 

w (ft ) = 1/4tt (26) 
and adopt Eq.(22) as the solution for the conditional probability. Hence Eq.(25) yields, 

m=—l I'm' 

x Jy (Qo) Yf* (Qo) dn J Yr (Q) Yf* (fi) dVt (27) 
6 



Using orthogonal property 

J Y™ (fi) Yf* (Q) dQ = 5 w 5 mml . (28) 
Hence, rotational correlation function (27) is written in terms of Mittag-Leffler function 

$ l (t)=$ l (Q)E a [-dJ(l + l) t a ] (29) 

where normalized factor of Eq.(29) is given as (0) = 4tt/(2/ + 1) . Eq.(29) states simply 
that the rotational correlation function, starting from the value unity at t — 0, decay non- 
exponentially in time with a relaxation time n that inversely proportional to rotational 
diffusion constant d a : 

Ta^idJil+lT 1 . (30) 

The rotational correlation function (29) is valid for arbitrary number I, and it has inter- 
esting property due to behavior of the Mittag-Leffler functions lies in the observation that 
is interpolates between an initial stretched exponential (i.e. KWW) behavior 



(t) ta exp 

and a long-time inverse power-law pattern 

$z (t) « 



r aiZ r (1 + a) 



(31) 



T a it 



(32) 



r(i-a) 

In conclusion, it is seen that relaxation functions Eq.(31) and (32) are compatible with 
Eqs.(2) and (3) for complex disordered systems, respectively. 

While drafting the problem we have disregarded the inertial of free motion governed by 
the rotational kinetic energy among the successive collision. Indeed, the rotational Brownian 
motion model for complex systems make sense only when friction r\ a is large i.e. the collisions 
are very rapid. There the rotational jumps were imagined to occur by large and arbitrary 
angles as opposed to the present instance in which only small- angle jumps are considered. 

In this study, we have analytically carried out the rotational relaxation function in terms 
of rotational correlation function for complex disordered systems based on rotational Brown- 
ian motion. To obtain rotational correlation function we have introduce generalized Fokker- 
Planck equations of fractional order, which generalizes the Stokes-Einstein-Smoluchowski 
relation, in consistency with the fluctuation-dissipation theorem. The introduction of the 



7 



Riemann-Liouville operator includes long-range memory effects which are typically found 
in complex systems, and consequently a single mode relaxes slowly in time, following the 
Mittag-Leffler decay. 

In conclude, we have shown that rotational Brownian motion in the complex systems such 
as spin glasses or dielectric materials, leads to KWW decays which indicates non-exponential 
relaxation. 
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